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Abstract 

New sufficient conditions for the characterization of dwell-times for linear impulsive systems are 
proposed and are shown to coincide with continuous decrease conditions of a certain class of looped- 
functionals, a recently introduced type of functionals suitable for the analysis of hybrid systems. This 
approach allows to consider Lyapunov functions that evolve non-monotonically along the flow of the 
system in a new way, enlarging then the admissible class of systems which may be analyzed. As a 
byproduct, due to the particular structure of the obtained conditions, the method is easily extendable to 
uncertain systems by exploiting some convexity properties. Several examples illustrate the approach. 
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I. Introduction 

Impulsive systems [QQ|, B2l, H, flU, 0, are an important class of hybrid systems in which 
the system trajectories admit discontinuities at certain time-instants. They arise in several fields 
like epidemiology Q, (H, sampled-data and networked control systems [|9]|, etc. Among the 
wide class of impulsive dynamical systems, we may distinguish systems whose impulse-times 
depend on the system state and those for which they are external and only time-dependent. 
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Linear systems of the latter class may be represented in the following form 

x{t) = Ax(t), t^t k , fee {1,2,...} 

x+(t) = Jx(t), t = t k , k G {1,2, ...} (1) 
x(to) = x 

where x,xo G W 1 are the system state and the initial condition respectively. The state x(t) is 
assumed to be left-continuous with x + (t) := lim s ^ t x(s) and the matrices A and J may be un- 
certain. The sequence {t\, t 2 , t 3 , . . .} is a strictly increasing sequence of impulse times in (t , oo) 
for some initial time to- The sequence is assumed to be finite, or infinite and unbounded. The 
distance between two consecutive impulse times is denoted by T k := t k+ i — t k with the additional 
assumption that e < T k , for some e > 0, which excludes the possible existence of a finite accu- 
mulation point. Given a set S C (0, oo), we define Is : = {{h, t 2 , . . .} : t k+ \ — t k G S, fceN} 
to characterize the sequence of impulse times. 

According to the matrices A and J, the system may exhibit different behaviors. In the case 
of periodic impulses with period T k = T, determining stability essentially reduces to study the 
Schurness of the matrix Je AT , which is a very simple problem. However, this formulation suffers 
from several critical drawbacks: 

1) the eigenvalue analysis is not extendable to the case of aperiodic impulses since the spectral 
radius is, in general, not submultiplicative; 

2) discrete-time Lyapunov approaches are thus needed, but they lead to robust LMIs with 
scalar uncertainties at the exponential, known to be complex numerically, but yet solvable 

OH; 

3) the extension to robust stability analysis is also difficult, again due to the exponential. There 
is no efficient solution, at that time, to handle matrix uncertainties at the exponential. 

The approach discussed in this paper aims at overcoming these problems. The method orig- 
inates from [fTTTl . [fT2l where an implicit but equivalent correspondence between discrete- and 
continuous-time domains is obtained by showing that discrete-time stability is equivalent to 
a very particular type of continuous-time stability proved using a new type of functionals, 
referred to as looped-functionals , lfT3~l . Since discrete-time stability focuses on the decrease 
of Lyapunov functions along the flow of a system at pointwise instants only, i.e. V(x(t)), 
t G {ti,t 2 )*3, • • •} U {oo}, its continuous-time extension V(x(t)), t G (0,oo) may then be 
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allowed to be nonmonotonic between consecutive pointwise instants. This hence authorizes the 
use of continuous-time Lyapunov functions that evolve non-monotonically along the flow of 
the system, enlarging then the class of admissible impulsive systems which can be possibly 
analyzed using existing approaches Q, fl4). Some results relying on a specific functional have 
been obtained for impulsive systems in O, |fl4l . Although leading to interesting results in some 
cases, selecting a functional however introduces some conservatism. 




Lyapunov Function V{x{t)) = x{t) T Px{t) 



Fig. 1. Continuous-time Lyapunov function V (plain) for Example |3.7| and the discrete-time envelopes (dashed). The discrete- 
time envelopes are monotonically decreasing while the continuous-time function is highly non-monotonic along the flow of the 
system. 

The current article develops a more generic framework based on an implicit characterization 
of a subset of the admissible looped-functionals. Since no particular functional is chosen in this 
paper, the obtained results are in essence more computationally-oriented than those in O, lfT4ll 
where functional are specifically chosen. The obtained conditions are expressed as robust convex 
optimization problems, conveniently depending on the inter-impulses interval distance and on the 
system matrices. They are also devoid of exponential terms and can be then efficiently extended 
to the uncertain case, usually not considered in the literature due to a lack of tractability. Thanks 
to the discrete-time nature of the stability notion, the approach is able to consider the jumps 
precisely by authorizing the continuous-time Lyapunov function to be strongly non-monotonic 
along the flow of the system; see Fig. [T] This feature allows to both characterize stability 
under periodic and aperiodic impulses. In the aperiodic case, stability under different dwell- 
times |fT51 , [fT6l , [0 are considered: minimal dwell-time T k E [T min , +oo), maximal dwell-time 
Tfc G (0,T max ) and ranged dwell-time Tk £ [T m i n ,T max \. The resulting feasibility problems are 
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then solved using polynomial programming, such as sum-of-squares [17] programming, very 
suitable for this kind of problems. 

The paper is structured as follows: Section |TT] is devoted to preliminary results on dwell-times 



for impulsive systems. In Section III alternative sufficient conditions for dwell-times are provided 
and applied in Section [IV] to uncertain systems. Finally, Section [V] is devoted to the presentation 
of Sum-of-Squares (SOS) programming used in the examples to enforce the conditions stated in 
the main results. Illustrative examples are treated in the related sections using SOSTOOLS lfT8ll 
and SeDuMi JT5I. 



The notation is quite standard. The sets of symmetric and positive definite matrices of di- 
mension n are denoted by § n and S2. respectively. Given two symmetric matrices A, B, A y B 
means that A — B is positive definite. For a square matrix A, the operator He(A) stands for the 
sum A + A T . The identity matrix of dimension n is denoted by I n . 

II. DWELL-TIME RESULTS 

A. Ranged Dw ell-Time 

The most general dwell-time notion is referred to as ranged dwell-time and is characterized 
in this section. Minimal and maximal dwell-times are corollaries of this result. 
Lemma 2.1 (Ranged Dwell-Time): Assume there exists a matrix P G §™ such that 

J T e ATe Pe M J -P^O (2) 

holds for all 9 G [T min ,T max ]. Then, for any impulse sequence in l[T min ,T max ], the system ([Tf is 
globally asymptotically stable. 

Proof: The condition ^ implies that V(x) = x T Px is a discrete-time Lyapunov function 
for the discrete-time system z(k + 1) = e ATk Jz(k), T k G [T min , T max ], since ^ is equivalent to 
the condition V(x(t k+1 )) - V(x(t k )) < for all t k+1 - t k G [T min ,T max ] and all x(t k ) G W 1 . 
This implies in turn that ||x(tfc)|| 2 — > as k — > oo. It remains to prove that the convergence of 
||#(£fc)||2 to implies the convergence of ||a;(£)||2 to when t — > oo. Defining t = t k + r, we 
have x(t) = e Ar Jx(t k ) and hence 

V(x(t)) < fi\\x{t k )\\l 

where fi = sup s6[Tm . n>Tmax] \ max (J T e AT T P e Ar J). Thus, when ||x(t fe )|| 2 -)■ then \\x(t)\\ 2 -)> 0. 
The proof is complete. ■ 
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The main difficulty of the above result lies in solving the robust semidefinite feasibility problem 
§2§ over the whole range of 9. The presence of exponential terms also prevents to easily derive 
a tractable generalization of this result to the uncertain case. 

B. Minimal Dwell-Time 

The minimal dwell-time result corresponds to the case when T min = T, for some T > 0, and 
Tmax — > oo. However, since the interval [T, +oo) is of infinite measure, it seems rather difficult 
to computationally check the condition ([2]) over such an interval. Fortunately, it is possible to 
provide a finite dimensional alternative result. 

Lemma 2.2 (Minimal Dwell-Time): Assume that for some given T > 0, there exists a matrix 
P e §™ such that 

A T P + PA -< (3) 



and 



jT e A T Tp e AT j _ p ^ o (4) 



hold. Then, for any impulse sequence in I[t,+oo> me system (|TJ) is globally asymptotically stable. 
Proof: The goal is to show that the conditions of Lemma 2.2 implies that J T e ATe Pe M J — 



P -< for all 6 e [T, +oo). A Taylor expansion of $(0) := e AT0 Pe Ae around 6 > yields 

$(# + S) = $(0 O ) + e ATf>0 {A T P + PA)e Ae °5 + o{6). 

Hence, condition ^ implies that the map is strictly decreasing, in the sense that $(# +e) -< 
®(9 ), for any 9 > and e > 0. Combining this with Q yields 

jT e A T (f + e) pe A(f + e)j_ p ^ Q 

for all e > 0. The proof of convergence of the state x(t) is finally proved as in the proof of 
Lemma 12.11 ■ 

C. Maximal Dwell-Time 

The maximal dwell-timeQ result can be obtained by letting T max = T, for some T > 0, and 
Tmin — > in Lemma [2T[ In this case, the interval (0,T] is of finite measure and the problem 

'This is also referred to as reverse dwell-time in (5) 
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of infinite measure does not hold. It is hence enough to consider the condition ([2]) over this 
interval. 

It is however possible, under the restriction that A is anti-Hurwitz, to reduce the condition to 
a finite-dimensional one, similarly as for the minimal dwell- time. 

Lemma 2.3 (Maximal Dwell-Time): Assume that for some given T > 0, there exists a matrix 
P G SI such that 

A T P + PA y (5) 

and 

jT e A T fp e Afj _ p ^ o (6) 

hold. Then, for any impulse sequence in I(o,t]> me system (|TJ) is globally asymptotically stable. 
Proof: The proof is similar to the one of Lemma |2.2[ ■ 

III. Alternative dwell-times characterization 

The main and common drawback of the results of Section [II] is the presence of exponential 
terms, preventing any extension to robust dwell-time characterization. Efficiently dealing with 
matrix uncertainties at the exponential is indeed an open and difficult problem. Checking the 
condition ([2]) over a finite-length interval is also computationally demanding, again due to the 
exponential terms. 

We provide here a way for approximating these LMI constraints involving exponential terms 
by exponential-less ones. These conditions are expressed as robust LMIs containing an addi- 
tional infinite-dimensional decision variable satisfying a certain boundary condition. It is also 
emphasized that the results can be interpreted as the continuous-time decrease condition of a 
very general looped-functional. 

A. Ranged Dwell-Time: 

As in the previous section, we start with the most general ranged dwell-time result: 
Theorem 3.1: The system ([[]) is globally asymptotically stable for any impulse sequence in 
\T min ,T max \ if there exist a matrix P e §'", a scalar e > and a matrix function Z : [0, T max ] x 
[T m i n , T max ] — >• §> 3r \ differentiable with respect to the first variable, verifying 

Y?Z(0, 9)Y 1 - Y 2 T Z(9, 9)Y 2 = (7) 
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for all 9 G [T min , T max ] where 





j 










In 


n = 


In 





, Y 2 = 












4 







4 



and such that the parameter dependent LMI 

tf(0) + He < 



holds for all r G [0, 9] and G [T min , T max ] with 

0(A T P + PA) 





.4 








T s 


< 











Z(t,6) > 
















d_ 



(8) 



¥(0) := 
















J T PJ - P + el n 



(9) 





Moreover, in such a case, the quadratic form V(x) = x T Px is a discrete-time Lyapunov 
function for the discrete-time formulation of system ([T]) and Lemma 2.1 holds with the same 
matrix P. 



Proof: To see the implication of Lemma 2.1 first pre- and post-multiply ([8]) by £fc(r) : = 
co\(x(t k + r),x(t k ),x(t k+1 )), where x(t k + r) = e Ar Jx(t k ), r G [0,T k ], t k+1 = t k + T k , to obtain 



Urf^(T k )Ur) + ^ [UrfZ(T,T k )Ur)} < 



(10) 



where we have set 9 = T k since we are considering the interval [t k ,t k+ i\. The second term of 



(10) has been obtained using the fact that 

d 



dr 



A 





(ID 



Integrating (10) from to T k we get 

% := J^[x(t k ) T (J T PJ~P + eI n )x(t k ) 



+ T k ^-V(x(t k + r)) 



dr 



(12) 



+Ck(T k ) T Z(T k ,T k )UT k ) - C(0) T Z(0,T fe )4 + (0) 
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where V(x) = x T Px. Noting that 



4 + (0) = Y 1 



x(t k ) 



x(t k ) 
x(t k+1 ) 



(13) 



we can see that the constraints make the two last terms of (12) equal to 0, and hence we 
have 



7] k = T k [V(x + (t k ))-V(x(t k ))+e\\x(t k )\ 
+V(x(t k+1 ))-V(x + (t k ))} 
= T k [V(x(t k+1 )) - V(x(t k )) + e | \x(t k ) 



(14) 



Since (10) is nonpositive, then r\ k < for all k and hence the feasibility of the conditions of 



Theorem 3.1 implies the feasibility of the conditions of Lemma 2.1 This concludes the proof. 



The above theorem provides then a sufficient condition for the asymptotic stability of linear 
impulsive systems for any arbitrary impulse sequence in ^-[T min ,T ma x]- The condition (|8j) can 
be understood as the continuous decrease over r G [0, 9} of the following family of looped- 
functionals parameterized by 9 E [T min ,T max ]: 



W := x{t k + r) T Px(t k + r) + &(r) 2 Z(r, 9)^ k {r) + rx(t k ) T (J T PJ -P + e I)x{t k ) (15) 



defined for r e [0,0] and satisfying the looping condition 

£+(0f Z(0,^+(0) = U0) T Z{9,9)U0)- 



(16) 



The obtained conditions however take the form of an infinite dimensional semi-infinite feasibility 
problem, which may be hard to solve. When the matrix function Z(t, T) is chosen in the set of 
matrix polynomials, sum-of-squares techniques can be applied in order to solve this problem. This 
framework allows to easily consider the equality constraint (looping-condition) as an equality 
constraint on the polynomial coefficients. One of the advantages of the proposed alternative 
formulation ([8]) is to provide conditions which are structurally well suited (convex) for the 
generalization to the uncertain case. Several dwell-time results for uncertain systems are derived 
in Section IW] 

Example 3.2: Let us consider the system ([[]) with matrices jH 

(17) 





-1 0.1 






1.2 





A = 


1.2 




J = 





0.5 
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Since both A and J have unstable eigenvalues, the system is unstable for sufficiently small and 
large periods, and cannot be studied using many of the existing approaches flU, [|4]|. In order to 



determine the lower and upper bounds for the inter-impulse intervals, Theorem 3.1 is applied and 
leads to estimates of Table [TJ It is important to stress that, in the periodic case, the minimal and 
maximal periods are given by 0.1824 and 0.5776, showing that the same bounds are achieved in 
the aperiodic case and that the proposed result is nonconservative for a multivariate polynomial 
Z(r,T) of degree 3. 



order of Z 


Tmin 


T-max 


1 


0.2040 


0.5672 


2 


0.1824 


0.5774 


3 


0.1824 


0.5776 


4 


0.1824 


0.5776 



Estimates of the admissible inter-impulse intervals for the aperiodic system of ExampleI3.2I 



B. Minimal Dwell-time 



The exponential-less approximation of the conditions of Lemma 2.2 is given below: 
Theorem 3.3 (Minimal Dwell-Time): Assume that for some T > 0, there exist a scalar e > 0, 

a matrix P 6 8" and a continuously differentiable matrix function Z : [0,T] — > § 3n satisfying 

YfZfflYi - Y?Z{f)Y 2 = and such that the LMIs 

A T P + PA -< 



(18) 

f(A T P + PA) ol 

J T PJ - P + el n 



hold for all r G [0,T]. Then, for any impulsive sequence in Ipf )+00 ), the system (111) is globally 







.4 








T x 


+ He < 













Z(r) > 


















d 



asymptotically stable and Lemma 2.2 holds with the same matrix P. 
Proof: The proof follows from Lemma [2T2] and Theorem 3.1 



Example 3.4: Let us consider the system ([T]) with matrices (61 



A 



1 



J 



2 1 
1 3 



(20) 
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Since A is Hurwitz, the minimal dwell-time concept applies and using Lemma |2.2| we get the 
value 1.14053 with the matrix 

' 0.6351 0.0601 
0.0601 0.3649 



P 



Considering then Theorem 3.3 with polynomial Z(r) of degree 3, we use the toolbox SOSTOOLS 
lfT8ll for sum-of-squares programming to enforce the conditions of Theorem |3.3| and we get the 



value 1.14054, which is very close to the one computed with Lemma 2.2 In such a case, the 
following matrix is obtained 

0.3443 -0.1102 
-0.1102 0.2285 

Estimates of the minimal-dwell-time for different degrees for Z are given in Table |nj 



order of Z 


Train 


1 


4.3255 


2 


1.1407 


3 


1.14054 


4 


1.14054 


5 


1.14054 



order of Z 




1 


36.3071 


2 


4.0705 


3 


2.6893 


4 


2.2966 


5 


2.2209 


6 


2.2010 


7 


2.1979 


8 


2.1974 


9 


2.1983 


10 


2.1990 


11 


infeasible 



TABLE II 

Estimates of the minimal dwell-time for Example 
\3A\ 



TABLE III 

Estimates of the minimal dwell-time for Example 
H31 



Example 3.5: Let us consider now the system ([T]) with matrices flU 



A 



100 
-1 



J 



-0.9 
0.9 



(21) 



In the T-periodic case, this system exhibits a quite complicated behavior for small T's, alternating 
stable and unstable regions. Since A is Hurwitz, we can determine the minimal dwell-time using 
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Lemma [272] and we get the minimum dwell- time value 2.1254. Theorem |3.3| leads on the other 



hand to the estimates of Table III We can see that, in this case, a gap still persists and the 
smallest one is obtained with a polynomial Z(r) of degree 8. Numerical problems arise when 
using polynomials of degree larger than 10, resulting then in infeasible problems. 

The above examples show some limitations of the approach when the convergence of the 
sequence of polynomials Z(r) of increasing orders, towards a suitable function, is quite poor. 
In such a case, a high order polynomial would be needed to obtain accurate results, but this 
would inexorably lead to feasibility problems involving a large number of decision variables 
and suffering from numerical problems. It is indeed well-known that SOS -programming scales 
very poorly with the problem complexity, and it is a matter of future research to see how the 
complexity of the problem addressed in this paper can be reduced. 



C. Maximal dwell-time 

A general maximal dwell-time result can be straightforwardly obtained be setting T n 







in Theorem 3.1 When the matrix A is anti-Hurwitz, we have the following simpler result, 



corresponding to Lemma 2.3 



Theorem 3.6 (Maximal Dwell-Time): Assume that for some T > 0, there exist a scalar e > 0, 
a matrix P £ and a continuously differentiable matrix function Z : [0,T] — >■ S> 3n satisfying 
Y^Z(0)Y 1 - Y?Z(f)Y 2 = and such that the LMIs 



A T P + PAyO 



f(A T P + PA) 

J T PJ - P + el n 









A 








T x 


+ He < 













Z(r) > 


















d 
dr 



(22) 



Z{t) d (23) 



hold for all r £ [0, T\. Then, for any impulsive sequence in I(o,t]> the system dlh is globally 



asymptotically stable and Lemma 2.3 holds with the same matrix P. 



Proof: The proof follows from Lemma [23] and Theorem 3.1 



Example 3.7: Let us consider the system ([T]) with matrices jH 



.4 



1 3 
-1 2 



J 



1/2 
1/2 



(24) 
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Since the matrix A is anti-Hurwitz and J is Schur, Theorem |3.6| may be applied and leads to 



the results of Table IV For comparison, Lemma 2.3 yields the value 0.4620, which is recovered 
by choosing Z(r) as a polynomial of order 3. 



order of Z 




1 


0.3999 


2 


0.4613 


3 


0.4620 


4 


0.4620 


TABLE IV 



order of Z 




1 


0.1067 


2 


0.1072 


3 


0.1072 



Estimates of the maximal dwell-time for Example 
IP 



TABLE V 

Estimates of the maximal dwell-time for Example 
I4~2l 



IV. Robust dwell-time analysis 

In this section, we extend the previous results to the case of uncertain systems with constant 
polytopic-type uncertainties, that is 

AeA:=co{A u ...,A Na }, Je J := co {J u ...,J Nj } (25) 

for some integers Na,Nj > 0. The system matrices can hence be alternatively written as 

Na Nj 

A = a A,iA, and J = ^ a J,i J i 

i=l i=l 

where both a a '■= col^cr^i) and <7j := colder belong to the unit simplex. The exponential 
terms in the usual dwell-time conditions have been major obstacles to the extension of the results 
to uncertain systems. Thanks to the proposed alternative formulations, considering uncertain sys- 
tems is henceforth possible. The theorem below provides a result concerning the characterization 
of ranged dwell-time for uncertain systems: 

Theorem 4.1: Assume there exist a matrix P G §™, a scalar e > 0, Nj matrix functions 
Zi : [0,T max ] x [T min ,T max ] -> § 3n , i = l,...,Nj, differentiable with respect to the first 
variable, verifying the condition (5^)^(0, 9)Yl - Y?Zi(6,6)Y 2 = for all 6 e [T min , T max ], 
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l,...,Nj where 





Ji 










In 




In 





, Y 2 = 


In 










In 







4 



and such that the parameter dependent LMIs 



9(A]P + PAj 













JfPJi- P + el n 






f 


A, 








T s 


+ He < 













Zi(r,9) > 


















d_ 



(26) 



N 



A- 



hold for all r £ [0, T], T £ [T min , T max ] and i = 1, . . . , Nj and j = 1, 

Then, the system ([T])-(25) is asymptotically stable for any impulse sequence in \T mm ,T max \ an d 
the following equivalent statements hold: 

1) The quadratic form V(x) = x T Px is a discrete-time Lyapunov function for the uncertain 
time- varying discrete-time system z(k + 1) = e ATk Jz(k) for any T k £ [T min ,T max ] and 



any (A, J) E Ax J. 
2) The LMI 



J T e AT9 Pe M J-P^0 



(27) 



holds for any 9 £ [T min , T max ] and any (A, J) £ Ax J. 



Proof: Multiplying (26) by <Jaj, summing over j = 1, . . . ,Na and integrating over [0,T] 



leads to the equivalent LMI Jj e ATT Pe AT .J { - P -< 0. Since the LMI holds for all z = 1, . . . , Nj, 
then using a Schur complement, it is possible to show that this is equivalent to J T e ATT Pe AT J — 
P -< 0, for all (A, J) £ Ax J. The proof is complete. 



Let us illustrate the above result with the following example: 
Example 4.2: Let us consider the polytopes 



A := co 



1 3 
-1 2 



3 1 
-2 4 



J := co 



1/4 
2/3 



2/3 
1/2 



(28) 



This uncertain system satisfies the maximal dwell-time conditions since the matrices in A are 



anti-Hurwitz. We apply Theorem 4.1 with Ziij) and adding the conditions AjP + PAj y 0, 
j — 1, . . . , Na- The computed dwell-time estimates for different degrees for the polynomial Zi 
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are gathered in Table |V} For comparison, the conditions of Lemma |2.3| are gridded and their 
solving yields the value 0.1072, which shows that the proposed approach gives the exact value 
for a polynomial of degree 2. Note, however, that solving the set of gridded LMIs is much longer 
than solving the proposed result. Using gridded LMIs is also less accurate since the LMIs are 
not checked over the whole set A. 

V. General presentation of Sum of Squares as Algorithmic Tool 
The methodology we use to implement the conditions of Theorems |3.1[ 3.3 3.6 and 4.1 



is based on the sum-of-squares (SOS) decomposition of positive polynomials. When applying 
this methodology, we assume that all matrix functions are polynomial, can be approximated by 
polynomials, or there is a change of coordinates that renders them polynomial. 

Denote by M[y] the ring of polynomials in y — (y%, . . . , y n ) with real coefficients. Denote by 
E s the cone of polynomials that admits a SOS decomposition, i.e., those p E for which 
there exist hi E M.[y],i = 1, . . . , M so that 

M 



p(y) = ^2^(y). 



i=i 

If p(y) E E s , then clearly p(y) > for all y. The converse is however not true, except in some 
very specific cases. The advantage of SOS is that the problem of testing whether a polynomial 
is SOS is equivalent to solving a semidefinite program [T71, hence worst-case polynomial-time 
verifiable. Note that on the other hand, checking whether a polynomial is positive is NP-hard. 
When extended to matrix polynomials, positivity is substituted by semidefinite positivity. The 
semidefinite programs related to SOS can be formulated efficiently and the solution can be 
retrieved using SOSTOOLS lfT8l . which interfaces with semidefinite solvers such as SeDuMi 

m. 

Consider now the conditions in Theorem 13.31 which take the form 



L(r)<0, teS:=[0,T], (29) 
where L(r) E § 3n and S is a semialgebraic set described by polynomial inequalities: 

S = {s E R : gi{s) > 0, % = 1, . . . , M}, 



where the g/s are polynomial functions. In order to test condition (29), we can apply Posi 



tivstellensatz results, such as the one in [|20l . which allow to test positivity on a semialgebraic 
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set using SOS programming. Specifically, Condition ( [29] ) holds if there exist SOS polynomials 
Pi(r, y), such that 



M 



i=l 

Intuitively, the above condition guarantees that when r G S, we have L(r) < — J2iL\ 9i( r )Pi( r i v) < 
since gi > and > 0, and therefore L(r) < for those r. 

VI. Conclusion 

Alternative sufficient conditions for the dwell-time characterization of linear impulsive systems 
have been developed. They can be interpreted as continuous-time decrease conditions of a certain 
class of looped-functionals. Ranged, minimal and maximal dwell-times have been considered. 
The obtained conditions are expressed as robust semidefinite programming problems which 
may be solved using polynomial techniques, such as sum-of-squares. Thanks to the convenient 
structure of the conditions, the results have been extended to uncertain systems. Several examples 
illustrate the efficiency of the approach. Future works will be devoted to the derivation of 
necessary and sufficient alternative conditions, stabilization criteria and the extension to nonlinear 
systems. 
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